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By proposing device-independent protocols, S. Pironio et al. [Natnre 464, 1021-1024 (2010)] and 
R. Colbeck et al. [Nature Physics 8, 450-453 (2012)] proved that new randomness can be generated 
by using perfectly free random sources or partially free ones as seed. Subsequently, Li et al. [Phys. 

Rev. A 84, 034301 (2011)] studied this topic in the framework of semi-device-independent and 
proved that new randomness can be obtained from perfectly free random sources. Here we discuss 
whether and how partially free random sources bring us new randomness in semi-device-independent 
scenario. We propose a semi-device-independent randomness expansion protocol with partially free 
random sources, and obtain the condition that the partially free random sources should satisfy to 
generate new randomness. In the process of analysis, we acquire a new 2-dimensional quantum 
witness. Furthermore, we get the analytic relationship between the generated randomness and the 
2-dimensional quantum witness violation. 

PACS numbers: 03.67.Ac, 05.40.-a 


I. INTRODUCTION 

Perfectly free random bits have both theoretical and 
practical significance. In the aspect of theory, perfectly 
free random bits are beneficial for the foundation of phys¬ 
ical theory to establish symmetries [1]. In practical ap¬ 
plications, perfectly free random bits could be used in 
many important fields, especially in cryptography. Al¬ 
most all the security of cryptographic protocols depends 
on perfectly free random bits. For example, in the well 
known BB84 protocol [2], the security will be seriously 
limited once an eavesdropper uses partially free random 
bits to replace the perfectly free ones [3]. 

Recently, the studies of device-independent (DI) and 
semi-device-independent (SDI) protocols have attracted 
a lot of attention. Here, DI means that no assumption 
is made on the devices used to perform protocols [4]. 
Subsequently, M. Pawlowski introduced the concept of 
SDI meaning that the devices in protocols are nonchar- 
acterized except the tight bound of the dimension of the 
potential required systems [5]. 

Randomness expansion is the protocol in which ran¬ 
dom sources are used as seed to produce new random¬ 
ness. Recently, R. Colbeck proposed a DI randomness 
expansion protocol based on the tripartite GHZ-type en¬ 
tangled states [6] and S. Pironio et al. proposed the 
protocol based on Bell inequality violation [7]. These 
results demonstrated that perfectly free sources can be 
expanded in the framework of DI. In 2012, R. Colbeck et 
al. showed that new randomness can also be obtained by 
using partially free bits as seed in the framework of DI 
(more precisely, the partially free bits can be amplified 
to make perfectly free ones and this process also is a DI 
randomness amplification protocol) [1]. Subsequently, Li 
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et al. studied this interesting topic in the framework of 
SDI and proved that new randomness can be produced 
from perfectly free sources by presenting SDI random¬ 
ness expansion protocols [8, 9]. Therefore, whether and 
how partially free sources bring us new randomness in 
the framework of SDI is a problem about which people 
may be curious. 

Here, we demonstrate that new randomness can be 
generated from partially free sources in the SDI scenario 
by proposing a SDI randomness expansion protocol with 
partially free sources. Different from the assumption that 
ei = £2 in the Ref. [1], we consider a more general 
case,where ei = £2 is not strictly required in our proto¬ 
col, and obtain the condition that £i ,£2 should fulfill to 
generate new randomness (the choices of states and mea¬ 
surements are derived from ei-free source and 62 -free 
source, respectively). A new 2-dimensional quantum wit¬ 
ness is gained in the process of randomness certification. 
Furthermore, the analytic relationship between the gen¬ 
erated randomness and the 2-dimensional quantum wit¬ 
ness violation is acquired. 

This paper is structured as follows. In Sec. H, we recall 
the definition of partially free random sources and intro¬ 
duce a SDI randomness expansion protocol with partially 
free sources. In Sec. HI, the condition which partially 
free sources should satisfy to generate new randomness, 
and certification parameters are obtained. In Sec. IV, 
the analytic relationship between the generated random¬ 
ness and 2-dimensional quantum witness violation is con¬ 
cluded. In Sec. V, we summarize our results. 

II. MODEL DESCRIPTION 

In order to better explain our theory, first of all, we 
give a detailed definition of partially free random sources 
in this section. 

Let V be a variable, considering its causal structure 
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FIG. 1. SDI randomness expansion with partially free sources. 
The dashed line represent that the hidden variable A may be 
correlated with these parts. Our protocol consists of two black 
box, which do not contain entanglement, in safe area. 


in relativistic space time, we call a variable A cannot be 
caused by X if A are not in the future lightcone of X. 
Denote the parameter A as the set of variables which 
cannot be caused by X and are interested in our devices. 
The variables in A may be provided by an eavesdropper 
or a higher theory [ 1 ]. 

Definition 1. A variable bit X is called e-free bit, 
e < 5 , if it satisfies |F’(0|A = A) — ^1 < e for all A S A. 
Particularly, X is called perfectly free bit as £ = 0. 

In this paper, we say that bits are picked according to 
e-free source if each bit is e-free and independents of 
other bits. 

Secondly, we introduce a SDI randomness expansion 
protocol with partially free sources based on n —>■ 1 quan¬ 
tum random access codes (QRACs) (see Fig 1). Based 
on the typical causal structure of our protocol [ 1 ], we as¬ 
sume that A may be correlated with two sources of weak 
randomness, the states prepared by Alice and the mea¬ 
surements performed by Bob. 

A detailed description of our scenario is described as 
follows: Alice picks n bits a = ooai • • • a„_i according to 
ei-free source Si and encoded to 1 qubit Pa,\, then Alice 
sends it to Bob via quantum channel. Bob performs two 
dimensional measurement = 0 , 1 } decided by 

?/ = 0 , 1 , • • •, n — 1 which is picked according to £ 2 -/ree 
source S 2 , and emits the measurement outcome b. In 
particular, there is not entanglement in the devices. 

In this paper, we construct the 2-dimensional quantum 
witness using the expected success probability which is 
different from the Ref. [11] and draws better conclusions. 

The expected success probability for the scenario is 

E = '^ P{a, y)P{b = |a, y) = '^ P^Ex, (1) 

a,y X 

where Ex = Y,a,y = ay|a,j/. A) and 

P{b\a,y,X) =tr{pa,x’M^^x)- 

Probability distribution of P{a,y,b) can be estimate 
by repeating the procedure many times, the value of E 
can then be estimated. 


Thirdly, we introduce the definition of the min-entropy 
function: 

Hoo{B\A, Y,A) = - log 2 max V P(A)P( 6 |a, y, A) (2) 

a.y.b.X 

aga 

to quantify the randomness of the measurement outcome 
for the scenario with the set A. 

Here the SDI randomness expansion with partially free 
sources based on 2 —1 QRAC is primarily discussed. 
The feasible region and the randomness certification of 
our protocol are explored in next section. 

III. FEASIBLE REGION AND RANDOMNESS 
CERTIFICATION 

In the DI randomness amplification proposed by R. 
Colbeck et al. [1], only one case of ei = £2 is discussed, 
and the relationship between quantum dimension witness 
and the min-entropy bound cannot be given as there are 
infinite parameters needed to be considered. In this sec¬ 
tion, we relax the assumption of £1 = £ 2 . Namely, the 
random resources of Alice is actually not required to be 
same as Bob’s, and obtain the feasible region. The good 
partially free sources are quite precious resource, our set¬ 
ting benefits to allocate partially free sources more rea¬ 
sonably and effectively. On the other aspect, the figure of 
the relationship between 2 -dimensional quantum witness 
violation and the min-entropy bound will be obtained 
through an optimization process. 

Definition 2. If there exists a protocol about SDI ran¬ 
domness expansion with partially free sources where Al¬ 
ice and Bob have the ei-free source Si,i = 0,1, respec¬ 
tively, and new randomness is certified, the pair (£i,£ 2 ) 
is called a feasible pair. The Feasible Region R of SDI 
randomness expansion with partially free sources is the 
set of all feasible pairs (£i,£ 2 ). 

It is evident for any A G A that the randomness ex¬ 
tracted from the outcome b will reduce to 0 with the 
increase of the distance between the probability distri¬ 
bution of P{a,y\X) and the uniform distribution on a,y. 
We assume that the eavesdropper attacks our devices in 
order to make our protocol get the least randomness and 
attempt not to led us finding that the random sources has 
been changed. To achieve his targets, the eavesdropper 
has to let 

\P{ai = 0|A) - = £i,i = 0,1. 

\P{y = 0\X)-^\ = e 2 , (3) 

for any A G A and 

P(a,y) = ^P(A)P(a,y|A) = i (4) 
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for any a £ {00,01,10,11}, y £ (0,1}. 

Without loss of generality, we can assume that there 
are only 8 hidden variables A^, fc = 0,1, • • •, 7 correspond¬ 
ing to 8 cases in Eq. (3). For the sake of convenience, 
let 

P(a, =0|Afc) = i + (-l)'=-ei,* = 0,l. 

P{v = m = \ + {-^f^e 2 . (5) 

where kokik 2 is the binary notation of k. It is easy to see 
that the eavesdropper can achieve Eq. (3) and Eq. (4) 
at the same time, see appendix A for the proof. 

Under the above attack of the eavesdropper, if for a 
pair (ei,S 2 ), 2-dimensional quantum witness violation 
still exist and the min-entropy is larger than 0 as 2- 
dimensional quantum witness violation reach its maxi¬ 
mum, then we can say that the pair {ei,£ 2 ) belongs to 
the feasible region R. 

Denote Ex^^c as the expected success probability with 
parameter Afe through a classical process. For any fc, 
the maximum value of can be obtained using the 

encoding map aoOi —>■ the bit decoding map 

0 —7> 0,1 —>■ 1 and the bit a(i_fc2) decoding map 0,1 —>■ 
and Ex^^c reach the same maximum value for 
any k. Obviously, the maximum value of E through a 
classical process is 

3 1 

Ec = — -l- 2^^! ^ 2 ) ~ ei£2- (6) 

Denote Exk,q as the expected success probability with 
parameter Afe through a quantum process, Ex^ appar¬ 
ently meet the linear relationship. For any k, to reach the 
maximum value of Ex^^q, as a general rule, every quan¬ 
tum state Pa,\k and positive operator valued measure 
(POVM) performed in the 2-dimensional 

space should be considered. Here, each mix state can 
be written as a convex combination of pure states. On 
the other hand, any POVM can be described as a con¬ 
vex combination of projective measurements, which in¬ 
clude projective measurements with rank 1, measure¬ 
ments {7,0} and {0,/} [12]. Different from the Ref. [8, 
9], we pinpoint that measurements {7,0} and {0,7} also 
need to be considered in our protocol. Nevertheless, we 
can prove that once measurements {7,0} or {0,7} are 
chosen, the relationship Ex^,q < Ec will be satisfied and 
is tight. In conclusion, here only pure states and projec¬ 
tive measurements with rank 1 need to considered. 

To visualize the pure states and projective measure¬ 
ments with rank 1, we consider Bloch sphere repre¬ 
sentation. Without loss of generality, set the Bloch 
sphere representation of measurement as 

-I’y.Afc} and Uo.Aj, = (1,0,0) for any k. Let the 
Bloch vector ra,Xk be the Bloch sphere representation of 
the pure state Pa,\k ■ 

For any pair (ei,£2), define 

t =_ 8£{(l + 4£i) _ 

1 + 16£t - 4£i - 64£t£i - ^ ^ 



FIG. 2. The feasible region R of SDI randomness expansion 
with partially free random sources is the yellow region under 
the green line but does not include the green line. Alice pick 
a according to £\-free source Si and Bob pick y according to 
£ 2 -free source S 2 , respectively. 

If t > 1, the optimal encoding-decoding strategy for 
any fc: Vi^Xk = (1,0,0) and r^.A^ = w^.aJ|| t’a.Aj|- 
Hence the maximum value of E is 

E = — + —£2 -h £1(1 — 2£2) < Ec- (9) 

If t < 1, the optimal encodin g-deco ding strategy for 
any k: Vi^Xk = ((-I)'=o+'=it, Vl - P,0) and ra,Xk = 
t’o.Afe/ll^^o.Afcll- After the simple analysis and calculation, 
Exk,q will reach the same maximum value and the max¬ 
imum value of E through a quantum process is 

^9 = 2 + 2 ^2 (10) 

For E = Eq, the min-entropy is 

where <5 = (I -I- 2 £2)/(I — 2£2). 

This implies that a pair (£i,£2) belongs to the fea¬ 
sible region R once it satisfies t < 1 , Ec < Eq and 
77oo(7?|A,y. A) > 0. Then the feasible region is ob¬ 
tained and demonstrated in Fig 2. Moreover, a new tight 
bound for 2-dimensional classical and quantum systems 
are given as Ec and Eq, respectively. Namely, a new 2- 
dimensional quantum witness is presented. See appendix 
B for a detailed calculation of Eq. (9)-(II). 

Next, for arbitrary pair (£i,£2) £ R, we begin to dis¬ 
cuss the min-entropy bound for a given expected success 
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probability E, which can be resolved by the following 
optimization problem: 

min H{B\A,Y,K) 

a,y,b,X 

7 

subject to : E = '^P{Xk)Ex^, 
k=0 

E\^ ^'^P{a,y\\k)P{h = ay\a,y,Xk), ( 12 ) 

a,V 

the optimization is carried out by quantum states 
Pa,A and POVMs {M^ chosen in 2-dimensional 

Hilbert space for a £ {00,01,10,11} and y £ {0,1}. 

After that we can estimate the min-entropy bound. 
Then true random numbers can be produced by a ran¬ 
domness extractor [13]. In fact, it plays an important 
role on many aspects that the min-entropy bound can 
be estimated as the analytic function of 2-dimensional 
quantum witness violation, such as security analysis of 
SDI randomness expansion [14]. 


IV. ANALYTIC FUNCTION 

In the SDI randomness expansion proposed by Li et 
al., the figure of the relationship between 2-dimensional 
quantum witness and the min-entropy is given, but the 
analytic relationship is not discussed. In this section, 
for arbitrary pair (ei,£ 2 ) G R, we explore the analytic 
relationship between 2-dimensional quantum witness 

7 

E = Y^P{Xk)E,^, (13) 

k=0 

and the min-entropy bound H{B\A,Y,A) = — log 2 P, 
where 

l/2+{t + d)/{2y^S‘^ + 2td+l) <p<l 

deduced from Eq. (11). 

We might take Aq as example. To depict a encoding¬ 
decoding strategy influenced by the hidden variable Aq, 
we extract two parameters {E\g,ma,Xa^y^b Pib\o,,y, Xq)), 
which can be regarded as points in the 2-dimensional 
coordinate system. 

For a given encoding-decoding strategy, E\g can be 
said as the convex combination of success conditional 
expected success probabilities obtained by pure states 
and projective measurements, and maxa,y,h P(61a, y, Aq) 
is not more than the convex combination of the maxi¬ 
mal guess probability obtained by the same pure states 
and projective measurements. That is, for a given value 
of E\g , the convex set composed of the realizable points 
achieved by pure states and projective measurements will 
provide a upper concave bound for maxa,y^b Pib\o,, V, Aq), 
denote the upper bound as pxg. We only discuss this 
upper concave bound in the following. 


Apparently, p\g can be viewed as a concave function 
of Exg. Denote pxg = C{Exg), C is a concave function. 
On the other hands, with the increase of Exg, the ran¬ 
domness generated by the quantum process will also be 
monotone increasing, as consequences C is a continuous 
and decreasing function. 

Fortunately, this discussion also applies to other hid¬ 
den variables Xk,k ^ 0. For any realizable point 
(Aao, maxa,y,h P(&la, y, Ao)), other hidden variables can 
realize through a single code. It is to say that other hid¬ 
den variables will reach the same bound px^ as pxg for 
Exk Exy: i.e., pxk ^(-^Aa,)- 

For the given E as indicated in Eq. (13), the lower 
bound of the min-entropy is 

7 

iL(H|A, Y,A) = -log 2 VP(Afc) max P{b\a,y,Xk) 

a,V,b,\k 

k—0 

7 

= -log2^P(Afe)Pv.. (14) 

fc =0 

Using the Jensen’s inequality, it is natural that if and 
only if Ex^. = Ex^, ,ky^k', the lower bound of min- 
entropy will be reached. Without loss of generality, we 
might take E = Exg , then p = pxg ■ The lower bound of 
min-entropy can be described as 

H{B\A,Y,K) = -\og^C{E). (15) 

The next work mainly describe the function C. Denote 
El as 

Ei=m&^{Exg-.C{Exg) = l}. (16) 

Obviously, we can deduce that C{E) = 1 as E < Ei 
according to the monotonicity of the function C. There¬ 
fore, the function C can be completely depicted once the 
one in the closed interval [Ei^Eq] is obtained and the 
value of El is determined, which only are discussed in 
the following. 

Based on the Refs. [8-10], the function C is monotropic 
as E £ [El , Eg]. Then the function C has a inverse func¬ 
tion denoted as C~^, i.e., E = C~^{p). The function 
C~^ can be obtained by the following optimization: 

max E = 'S^P{a,y\Xo)P{b = ay\a,y,Xo) 

a,y,b,\o ^' 

a,V 

subject to : max P(5ja, y, Aq) = p, (17) 

a,y,b,\o 

where 1/2 -I- (t -I- 6)/(2x/S^ + 2t6 + 1) < p < 1, the op¬ 
timization is carried out pure quantum states pa,x and 
projective measurements chosen in 2-dimensional Hilbert 
space for a G {00, 01,10,11}, y G {0,1} and A G A. 

Firstly, we focus on the value of E achieved by arbi¬ 
trary pure states and projective measurements with rank 
1 in the optimization (17). Fortunately, E can be viewed 
as a function G(ei, £ 2 ,p) determined by ei, £2 and p (see 
the appendix C). 
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Secondly, consider the optimization (17) achieved by 
arbitrary pure states and measurements {/, 0 } , { 0 , /} in 
optimization (17), only E = Ec ss p = 1 is obtained. 

If > Ec is established. The strategy with mea¬ 

surements {/, 0 } or {0,7} can be simulated by one with 
pure states and projective measurements with rank 1 , 
where = G'(ei,£ 2 , !)■ Here E = G(ei,e 2 ,p)- 

On the contrary, if < ^c, then we have to 

discuss the convex set which is composed of points 
(G(ei,e 2 , 7 ')) 7 ') and {Ec,l) to obtain the upper bound 
of E denoted as 7'(ei, £ 2,^)5 which is a function of £i ,£2 
and p. In fact, the points {E(ei, £ 2 ,p),p) also provide 
a lower bound of min-entropy. However, whether the 
bound is tight or not cannot be determined. 

Based on the above analysis, we have Ei = 
maxlTic, and for Ei < E < Eg, 

P ; ^ 21 p) ; '^f ^51,52 — 

^^^“\F(ei,£ 2 ,p) ifE,„,,<E,, 

where 1/2 -|- (t -|- 5)/{2\/5'^ + 2t6 + 1 ) < p < 1. 

In fact, the function G and E can describe the re¬ 
lationship between the min-entropy bound and the 2 - 
dimensional quantum witness violation in detail. 

Denote /3 = arccos(2p — 1), we have 

G(ei,e 2 ,p)= ^ max {Gi(ei, £ 2 ,p, a)}. 

[ 0,77 — 4 / 3 ] , 2 = 1,2 

The analytic functions Gi, G 2 are 

Gi(ei,£ 2 ,p,a) = ^ + i(i - - e 2 )[( 5 cos /3 

-f cos(/3 -f a) -7 /(£i, £ 2 ,P, a)]- (18) 

G2(£i,£2,P,a) = ^ + ^(^ - - £2)[5crcos/3 

-f cr cos(/3 -7 a) -f 5 f(£i, £ 2 , p, a)], 

where cr = (1 -I- 2£i)/(l — 2£i) {f,g see Eq. (C12), Eq. 
(C13), respectively). 

The function E can be depicted by the function G and 

F{£i,£2,p) = (G(£i,£2,Po) “ Ej (1 - p)/(I - Po) 

+ E, (19) 

as p > Po, 7 ^(£i,£ 2 ,p) = G(£i,£ 2 ,p) as p < po, where po 
satisfies 

G{£i,£2,Po) - Ec . G(£i,£2,p) - Ec 

- 1 -= min{-^-}. 

Po - 1 p p - 1 

In particular, for the cases of £i = £2 < 0.1358; £1 < 
0.2203, £2 = 0 and £1 = 0, £2 < 0.5, we have 

Esi,s2 ~ aiaxQ,^[o^7r—4/3 ]{I^i(£ij£ 2,15 *a)}- 

In the Ref. [ 8 ], the analytic relationship between the min- 
entropy bound and the 2 -dimensional quantum witness 
violation is E = G(0, 0, 

Many works only consider projective measurements 
with rank 1 just because they happen to satisfy > 
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FIG. 3. The relationship between the min-entropy bound and 
the 2-dimensional quantum witness for 0 < e < 0.1358, where 
the choices of states and measurements are derived from e- 
free sources Si and S 2 , respectively. 

Ec, such as the Refs. [ 8 , 9]. We might take £1 = £2 = £ 
as example. In fact, Eg,^ > Ec as £ < 0.12348. But 
Eej_c 2 < Ec as 0.12348 < £ < 0.1358. This shows that 
the situation of < Ec will occur with the increase 

of £ 1 , £2 and the measurements {/, 0 } and {0.7} must be 
taken into consideration. Furthermore, the relationship 
between the min-entropy bound and the 2 -dimensional 
quantum witness for £1 = £2 is demonstrated as the Fig.3. 

V. CONCLUSION 

We proved that partially free sources can bring us 
new randomness and proposed a SDI randomness expan¬ 
sion protocol with partially free sources based on 2 —>• 1 
QRAC. In our protocol, the condition that the partially 
free sources should satisfy to generate new randomness 
was gained without strictly requiring £1 = £2 (the choices 
of states and measurements are derived from ei-free 
source and £ 2 -free source, respectively). Furthermore, 
a new 2 -dimensional quantum witness and the analytic 
relationship between the generated randomness and the 
2 -dimensional quantum witness violation were obtained. 
In addition, the advantage of no containing entanglement 
which is introduced in the Ref. [ 8 , 9] also apply to our 
protocol. We conjecture that it can get better results in 
the SDI randomness expansion with partially free sources 
based on n —>■ 1 QRACs for n > 3. 
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Appendix A 

The assumption in Eq. (5) apparently satisfy Eq. (3). 
In addition, Eq. (4) is equivalent to 

P{a, = 0) = ^ P{Xk)P{ao = 0|Afe) = -,i = 0,1. 

k^O 

^ 1 

p(y = 0) = ^P(Afc)P(j/ = 0|Afc) = -. (Al) 

k^O 

The Eq. (Al) can be written as 

^ P(Afc)- E PiXk) = 0. 

fe=0,l,2,3 fc=4,5,6,7 

^ P(Afc)- E P(>^k) = 0. 

fc=0,1.4,5 fc=2,3.6,7 

^ P(Afc)- ^ P{Xk)=0. (A2) 

fc=0.2,4,6 fc=l,3.5,7 

They are the linear equations of Afc, A: = 0,1, • • •, 7. It is 

easy to see there must be many solutions to Eq. (A2) as 
there are three equations but eight variables. 


Appendix B 

The expected success probability with variable Xk is 
= ^P(a,y|Afc)P(6 = ay\a,y,Xk) (Bl) 

a,V 

for fc = 0,1, • • •, 7. 

We might take the Aq as example. Set the Bloch sphere 
representation of the pure state Pa,\o ■, projective measure 
Bloch vector ra,Ao, 

for a G {00, 01,10,11} and y G {0,1}, respectively. With¬ 
out loss of generality, let vq.Ao = (IjOiO). By the Ref. 
[10], we can know 

P{h\a, y, Ao) = tr{paM^l,\o) 

= ^(1 + r’a.Ao • (B2) 

where “ • ” denotes the inner product. 

With a small amount of calculation, we get 

+ 2/l^o)T’a.Ao • 

= ^ + \Y1 -f’(a|^o)ra,Ao • Va,Xo 
a 

<l + lT.P(^\^0)\\Va,xJ 

a 

(B3) 


where Wo,Ao defined as Eq. (8). Only the case of 
II ^o,Ao II7^ 0 is discussed in the following. If and only 
if t’o.Aq = I’a.Ao/ll^’o.Aoll, the Eq. (B3) can achieve the 
maximum value. 

Furthermore, set 9 is the angle between 4 ^o,Ao and 
WpAo. then 

II^OO.Aoll^ + ||l’01,Aoll^ = 1-1- 4:62- 
P(00|Ao) + P(Il|Ao) = i+2£^ 

P(01|Ao) + P(10|Ao) = i-2£2, (B4) 

where Alice and Bob have the Si-free source S'!, i = 0,1 
to pick a, y, respectively. 

With the knowledge of Eq. (B4), we have 

^P(a|Ao)||^a.Aoll = (i+2£?)||^O0.Ao|| 

a 

+ (2 “ 2ei)||i;oi,Aoll 

“ \J 2~^ + ^^2- (B5) 

If and only if || vqo.Ao ll/ll^’oi.Ao II = (1 + 4ef)/(l - 4ef), 
that is, the angle 9 must satisfies cos0 = t, t is defined as 
Eq. (7), the Eq. (B5) can achieve the maximum value. 

If t < 1 for a pair (ei,£ 2 ). Fortunately, let 9 = arccost, 
we will reach the maximum value of E\g denoted as 


rpmax 

P>'o,q 


1 

2 


1 

2 



-I- Sef + 252 + 32£{£2- 


(B6) 


At the same time, the maximum success probability with 
variable Aq can be obtained: 


P(5 = 0|00,0,Ao) 


^ f1 4- t + d 

2 ^ '^^/WTWTf 


(B7) 


If t > 1 for a pair (£i,£ 2 ), we have to set cos0 = 1, 
i.e., d = 0 to obtain the maximum value 


EZ%^ = ^ + ^52+5!(1 - 252). (B8) 


Appendix C 

For the sake of simplicity, suppose P = Pao p = 
P\o- 

To reach the maximum value of E satisfying a 
condition that the probability distribution satisfies 
maxa,y,f, P{h\a, y, Aq) = p = ^(1 -I- cos/3), we consider the 
question that which one of 16 probability P{b\a,y, Xq) 
should reach p, where 


{t + S)/\/ (5^ -|- 2t6 -I-1 < cos /3 < 1. (Cl) 

Apparently, it must be a guessing success probability 
and need to satisfy the following conditions: 
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(i) It must be the probability P(h = ao|a,y = 0,Ao) 
for a G {00,01,10,11}. Since P{h = ao|a,?/ = 0,Ao) > 
P{b = ai|a, y = l,Ao) can achieve a larger value of E 
concluded from P{y = 0|Ao) > P{y = l|Ao). 

{ii) The bolch vectors Va^Xg, Vy^Xg for all a,y are in 
a plane, ra Xg fall on the area between (—l)““i;o Aq and 

Then only four cases of P{b = ao|a,y = 0,Ao) = p 
for aG{ 00,01,10,11} need to be discussed on the strict 
precondition (ii). It is noteworthy that 


Assume Ooo.Ao = ^’oo,Ao/I|i’oo.AoII for a G [ 61 , 62 ]- In 
the case of a G [0, 61 ) 11 ( 62 , tt— 4/3], let /3 the angle between 
aoo.Ao and uq.Aq eis (3, where 61 = arcsin((5sin/3) —/3, 62 = 
min{ 7 r — 4/3, tt — arcsin((5 sin /3) — /3}. 

It’s worth noting that 61 < 02 for any pair of the fea¬ 
sible region. 

Based on the analysis of the above, we can get the 
analytic function of E in the case 1 and denote it as 
Gi(ei,e 2 ,P,a) shown as Eq. (20). With some calcula¬ 
tion, we can obtain 


P{b = ai|a,y = 1, Ao) < P{b = ao\a,y = 0,Ao) (C2) 

always be established no matter which one of P{b = 
ao|a, y = 0, Aq) reach p. 

Case 1: Let 

P(6= l|ll,y = 0,Ao) =p, (C3) 

that is, the angle between — Vq.Ao Etnd is P, suppose 

the angle between rn Aq and —Vi Xg is P+a. Since P(6 = 
l|ll,y = l,Ao ) < p, let a > 0 . Fortunately, for any 
a > 0, there is at least one choice can satisfy P{b = 
OjOO, y, Ao) < p which let roo,Ao along the same direction 
as VQ^Xg + '*^l,Ao • 

In order to ensure 

P(b = 0|01,y = 0, Ao) < p, 

P(6 = l|10,y = 0,Ao)<p, (C4) 

then Of < TT — 4/3 concluded from [tt — (2/3 -|- a)]/2 > /3. 
We confirm that the range of a is [0, tt — 4/3]. 

If the value of a is determinated, the angle between 
Vq.Ao and Vi,Ao is determined, then the vector fi^Ao is 
determined since we have set i^o.Aq = (!> 0, 0). 

Next only need to consider how to place Va^Xg for the 
purposes of reaching the largest E. 

Firstly, roo.Ao has been determined. 

Secondly, denote p as the angle between toi^Aq and 
^’o,Ao • p < 7r/2 is obtained derived from the Eq. (C2). 
To get a lager value of E, we want to set 

r’oi.Ao = ■*^01.Ao/ll«01Aoll, (C5) 


Gi(£i,e2,p, a) = ^P(a,y|Ao)P(6 = aj^|a,y, Aq) 

= ^ + \{\-£i?{\-£2)[5cofil3 

-hcos(/3-ha) -h/(£i,e 2 ,p, a)] (C7) 

where f(ei,S 2 ,p,a) is displayed as Eq. (C12). 

Case 2: Let 

P(6 = 0|01,y = 0,Ao) =p, 

that is, the angle between —i»o,Ao and roi,Ao is /3, set the 
angle between rgi^Xg and — Vi^Ao is /3 + a. Obviously, the 
range of a remains [ 0 , 7 r — 4/3]. By a similar way, denote 
E as G 2 (£i,£ 2 ,P, a) in this case, we have 

G 2 (si,e 2 ,p,a) = ^ + i(i -£i)^(i -£ 2 )[ 6 crcos /3 

-I- cr cos(/3 -I- a) -I- y(£i,£ 2 ,p, a)](C 8 ) 

The detailed description of g is displayed in Eq. (C13). 
Case 3: Let 

P(6=l|10,y = 0,Ao)=p. 

The angle between —vq^Xo and Tio.Aq is /3, set the angle 
between toi^Aq and —Vi,Ao is l3+a. The analytic function 
of E is equal to G 2 (£i, £2,P, a) concluded from P(a = 
10 ) = P(a = 01). 

Case 4: Let 


but must guarantee p > /3 for the sake of that Eq. (C4) 
will not be established, i.e., if Eq. (C5) want to be es¬ 
tablished, we must guarantee that 


tanp 


sin(2/3 -I- a) 

5 — cos(2/3 -I- a) 


> tan/3. 


(C 6 ) 


P (6 = 0|00,y = 0,Ao) =p. 

Set the angle between roopo and Vipo is /3-|-a. Denote 
E as G 3 (£i,£ 2 ,p, a) in this case and a G [ 0 , 7 r — 4/3]. By 
the same analysis as the case 1 , we get 


Using the knowledge of trigonometric functions, the 
Eq. (C5) is equivalent to sin(3/3-|-Q;) > 6 sin /3. Combining 
the condition a G [ 0 , 7 r — 4/3] and sin(3/3 -I- a) > 6 sin (3 
yields that the Eq. (C 6 ) can be established only for a G 
[ 01 , 02 ], where oi = maxjO,arcsin(i5sin/3) — 3 / 3},02 = 
TT — 3/3 — arcsin((5 sin/3). 

If a G [0, oi) U ( 02 , TT — 4/3], we have to set p = /3 to get 
a lager value of E. 

At last, by the similar way, suppose rio,Ao = 
■*^io,Ao/ll’^ioAol] fora G [ 01 , 02 ]. IfaG [0, oi)U(o 2 , 7 r- 4 / 3 ], 
suppose the angle between rio.Ao and — Uq.Ao is /3, 


G3(£i,£2,p,a) = Gi(£i,£2,p,a) (C9) 

for a G [0, 61 ] U [ 62 , TT — 4/3]. 

For a G [ 61 , 62 ], we have 

G3(£i,£2,P,a) = i -h i(i - - £2)[o-^(5cos/3 

-I- x/5'^ + 1 + 25 cos(2/3 -I- a) + 
cos(/3 -|- a) -I- A:(£i, £ 2 ,p, a)] 

where 



f{£i,£2,P,o:) = < 


(C12) 


' {26a + a'^6) cos /3 + cr^ cos(^ + a) — 2cr cos(3^ + a) 

(7^(5 cos /3 + (T^ cos(/3 + a) + 2a \/ 6'^ + 1 — 25 cos(2/3 + a) 
a^\/6'^ + l + 26 cos(2/3 + a) + 2a \/ 5“^ + 1 — 26 cos(2/3 + a) 
, a^ \/6^ + 1 + 26 cos(2/3 + a) + 25a cos jS — 2a cos(3^ + a) 


if a € [ 0 , ai) U {62, tt 
if a€ [ai,bi), 
if a € [61,02), 
if a € [02,62]. 


4 / 3 ], 


g{ex+2,P,a) 


' 5a cos P + a cos(/3 + a) + (<7^6 + 6) cos /3 — (a^ + 1) cos(3/3 + a) 

6a cos P + a cos(/3 + a) + (cr^ + 1) + 1 — 26 cos(2/3 + a) 

O'a/ 6^ + 1 + 26 cos(2/3 + a) + (a^ + l)-\/6^ + 1 — 26cos(2/3 + a) 

, O'a/ 6^ + 1 + 26 cos(2/3 + a) + (a^6 + 6) cos/3 — (a^ + 1) cos(3/3 + a) 


if a G [ 0 ,oi)U (62,77-4/3], 
if aG [oi,6i), 
if a G [61,02), 
if a G [02,62]. 

(C13) 


k{ei,e2,P,a) = 

I 2 a 62 + 1-26 cos(2,8 + a) //a €[61,02) 1Q2.O) 

\ 26 a cos P — 2 acos{ 3 P + a) if a G [02, 62]. 


concluded from Eq. (C9) and Eq. (Cll), where Gi and 
G 2 are shown as Eq. (C7) and Eq. (C 8 ). 


Therefore, for a G [ 61 , 62 ], 

Gi(£i, e2,p, a) > G3(£i,£ 2,P,a) (Cll) 
always is established. We have 

G(£i,£2,_p)= max {Gi(£i, £2,p, a)}. 

aG [0,77—4^] = 1,2 


In particular, arcsin (6 sin /3) < 0 and tt — 

arcsin (6 sin ,5) — P > n — AP can be obtained by the tool 
of Matlab as £1 = £2 < 0.1358. Then oi = 0,62 = ir — AP, 
and the set of [0, oi) U ( 62,77 — 4/3] does not exist. The 
calculation process will becomes much more simple. 
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